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ABSTRACT 

Proteins play an important role in cell functions. Protein structures are useful for determining its functions. 
Dihedral angle places an important role in determining the structure of proteins. The aim of encoding the angle of 
dihedral or torsion angle is done with respect to the plane of two angles. This paper explains encoding the dihedral angle 
of amino acids in a protein by applying the mathematical representation of determining dihedral or torsion angle with 
respect to the plane of two angles. By calculating the cross product and dot product of two vectors of amino acid the 
dihedral angle are determined. 

KEYWORDS: Protein Structures , Dihedral angles , Torsion Angles & Plane of Angles, PDB 



TRANS 

STELLAR 

• Joumal Publications • Research Consultancy 


Receiyed: Jan 10, 2017; Accepted: Feb 17, 2017; Published: Feb 22, 2017; Paper Id.: I JMPERD APR20 174 

INTRODU CTION 


“A dihedral angle is defined as the angle between two intersecting planes. In Proteins, the dihedral angles 
are formed between the planes of two sets of three backbone atoms which can \(/, § or % l , x 2 , x 3 etc.” The dihedral 
angle is used to represent the accurate conformation of protein structures along with the force field. The \\t and cp 
angles contribute to predicting the structure of protein. [1] Dihedral angles are mainly useful in predicting the 
mimetic of proteins in drug discovery. [2] 

METHOD 


To encode the dihedral angle of a protein, it’s important to represent the protein in a proper manner. 
The representation of amino acids is extracted from a PDB (Protein Data Bank) data file. Figure 1 shows PDB file 
which contains the information with PDB number, Protein name, Residue Number, X, Y and Z co-ordinates in 
space. 


Atom 

PDBNo 


Protein Residye No 

X 

Y 

Z 

Atom 

ATOM 

1 

N 

HIS A 

1 

49.668 

24.248 

10.436 

N 

ATOM 

2 

CA 

HIS A 

1 

50.197 

25.578 

10.784 

C 

ATOM 

3 

C 

HIS A 

1 

49.169 

26.701 

10.917 

C 

ATOM 

4 

0 

HIS A 

1 

48.241 

26.524 

11.749 

0 

ATOM 

5 

CB 

HIS A 

1 

51.312 

26.048 

9.843 

c 

ATOM 

6 

CG 

HIS A 

1 

50.958 

26.068 

8.340 

c 


Figure 1: PDB - File Information of Amino Acids with Each Atom 


These PDB files do not contain the Hydrogen information; hence Hydrogen is added to the representation 
after accessing each residue. After Hydrogen is added, each atom is given an Atom Number and stored in a matrix 
along with its co-ordinates information. Figure 2 gives the representation of 3 amino acids after adding Hydrogen 
atom. [3] 
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Three Amino Acids HIS, SER, GLN are taken and represented as illustrated in the Figure 2. Starting from the first 
Atom Nitrogen (N) which is assigned 0, Ca assigned atom number 1, etc. All the backbone of an AA’s are assigned the 
atom number, and then followed by the side chain molecules; last, the Hydrogen atoms are assigned Atom numbers. 
Similarly, all AA’s are assigned the atom number. 

• Representation of Dihedral Angles 

After all AA’s are represented with atom number, a need to represent the dihedral angle is done. Figure 3 shows 
the representation of dihedral angle. By considering the dihedral angle to be encoded to be present in the i th AA, its 
previous AA is taken as i-1 and next AA is taken as i+1. Cai has § and \|/ angles respectively. The plane of two angles are 
given as Ci_i — Nj — C ai and Nj — C ai — Q. 



Figure 2: Representation of 3 Amino Acids of a Protein 


Ci-i — Ni — Cai — Ci 
M?: Ni — Co.i — Ci — Ni+i 
X i:Ni — Cai— Cpi— Cyi 

Figure 3: Representation of Dihedral Angle 

• Encoding the Dihedral Angles 

To encode the dihedral angle, select the atom number as illustrated in Figure 4. In the following case by 
considering the Figure 2, atom number 19 is selected and corresponding co-ordinates from the PDB file are accessed as 
shown in Table 1. 


4 v 

<|>: O — Ni&-£ Caie 0 - C20 

Figure 4: To Encode the Dihedral Angle of Cal9 Atom 
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Table 1: Co-Ordinates of 4 Atoms from PDB 


Atom 

Amino Acid 

X 

Y 

Z 

C 2 

HISTIDINE 

49.169 

26.701 

10.917 

N lg 

SERINE 

49.788 

27.850 

10.784 

C a 19 

SERINE 

49.138 

29.147 

10.620 

C 20 

SERINE 

47.718 

29.251 

10.110 


• Mathematical Representation of Dihedral Angle 

Dihedral angles are the angles between two planes as shown in Figure 5, mathematically it can be calculated using 
the Cross product and Dot product of Yectors. [4] 



Figure 5: Representation of Dihedral Angles between Two Planes 

The Eq. used to calculate the vector n^and njis given below. 

. AB*BC 

m = , — >n^i 

1 |ab||bc| 

. . BC*CD 

Eq. (1) and (2) represents the vector AB and BC and the magnitude of vector AB and BC respectively. 
cos(7i — x) = -^2 

cos x — — nj -^2 


( 1 ) 

( 2 ) 

( 3 ) 

(4) 


The dihedral angles of the atom number 19 in Figure 4 can be represented using the two planes as Plane 1 with C2 
- N18 - Cal9 and Plane 2 as N18 - Cal9- C20. Using the Eq. (1) and Eq. (2), the dihedral angle can be calculated as 
shown in the Eq. (5). The calculation of vector and the distance is shown in Table 2 and illustrated in Figure 6. 
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Figure 6: Representing Dihedral Angle Ca for Amino Acid SERINE in the Given Protein 
• Dihedral Angle Calculation 

The Eq. 1 and Eq. 2 to fit the amino acid is represented in Eq. 5 and Eq. 6 respectively. 

n _ C 2 N 18 *Ca 19 N 18 

1 _ |CiNlS| |Ca 19 N 18 | 

_ Ca 19 N 18 *C 20 Ca 19 

2 |Ca 19 N 18 | |C 2 oCa 19 | 

Calculating the X, Y, Z and Distance is illustrated below, for the vector, and the Eq. are as follows. 
The corresponding values after calculation are given in Table 2. 

X = Ca 19 X — N 18 X 
Y = Ca 19 Y — N 18 Y 
Z = Ca 19 Z — N 18 Z 

Dist3Lnce,\C al9 N 18 \ is calculated as shown in the Eq. 10 


(5) 

( 6 ) 

(7) 

( 8 ) 
(9) 


I C a i 9 Ni 8 1 X * \C al9 N 18 \ X -\- \C al9 N 18 \Y * \C al9 N 18 \Y -\- \C al9 N 18 \ Z * \C al9 N 18 \Z 


( 10 ) 


Table 2: Calculation of Yectors and its Distances 


Vector 

X 

Y 

Z 

Distance 

^al9^18 

-0.650 

1.297 

-0.164 

1.460 

c 2 n 18 

-0.619 

-1.149 

0.133 

1.311 

^20^al9 

-1.425 

0.104 

-0.51 

1.317 

^al9^18 * C 2 N 18 

-2.040 

-1.009 

0.101 

2.279 

^al9^18 * C 2 N 18 

* ^a!9^18 

-0.015 

0.187 

1.549 

1.561 


• Calculating Cross Product and Dot Product of Dihedral Angles 


The Cross Product ofC al9 iV 18 * C 2 N 18 is calculated using the equation as follows and the result is placed in 

Table 2. 


X — (Qi9^i8 Y * C 2 N± 8 z) (C al9 N 18 z * C 2 N 18 y ) 


( 11 ) 
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Y — (Ocl9^18 z * ^2^18 x ) (0*19^18 x * ^2^18 z ) (12) 

z = ( 0*19^18 x * ^2^18 J 7 ) — (fial9^18 y * ^2^18 x ) (13) 

The Dot Product for n^ and nj is calculated using the Eq.5 and Eq. 6 is applied for the amino acid atoms as 

follows. 


nl 


C 20 C al9-( C al9 N 18* C2 N 18)* C al9 N 18 

||Q*19Yl8* C 2 N 18 \.\C al9 N 18 \\ 


n2 = 


C 20 C al9-{ C al9 N 18* C 2 N 18 ) 

l| C al9^18|-l C 2^18ll 


Finally Dihedral angle \|/ is calculated by 

n i V (nT*nl ) 

C05 lp = ===== , = 

V ((n2)*(n2) )+ (n2) * (n2 ) 

The Dihedral angle \|/ in the Eq. 16 can be rewritten as 


ip = Cos — 1 


V(nr*nl ) 


V ((n2)*(n2) )+ (n2) * (n2 ) 


¥ = 


\|/ * 180 


(14) 

05) 

06) 

(17) 

(18) 


“The dihedral angle determines the conformations around bonds which are rotating. The change in the dihedral 
angle is possible with only the distance between the lst and 4th atoms of the amino acids; the other 2 atoms would 
contribute only for chemical bond lengths and bond angles. Its value range from -180 to +180 degrees. The dihedral angle 
is considered to be +ve if a clockwise rotation is performed with the molecule and it will be -ve when an anti-clockwise 
rotation is performed with the molecule in its plane.” [5,6, 7] 


If (nl < 0), v|/ = 180 — v|/ else v|/ = — vj/ 


(19) 


RESULTS AND DISCUSSIONS 

The dihedral angle for C al9 N 18 is calculated by accepting X of C a i 9 as 49.138 and N 18 of Eq. 7, X as 49.788, 
resulting with -0.65, the same is displayed in the Table 2. Similarly same value of X, Eq. 8, Y and Eq. 9, Z is calculated. 
The distance Eq. 10,|C al9 N 18 | = 1.460. The cross product of C al9 N 18 * C 2 N 18 is calculated by accepting the values from 
Table 2 and Eq. 11, X is calculated as (1.297 * 0.133) - (-0.164 * -1.149) resulting with -2.040. Similarly same value of X, 
Eq. 12, Y and Eq. 13, Z is calculated. The dot product Eq. 14, nlis calculated by accepting the values of X, Y and Z as 
{(-1.425 * -0.015) + (0.104 * 0.187) + (-0.51 * 1.549)} / (1.561), resulting with the value 1.3156. Similarly, Eq. 15,n2is 
calculated resulting with value -0.508. Finally Eq. 16, cos vj/ = 0.9327 and Eq. 17 vj/ is the cos-1 with value 0.3689. This Eq. 
18, v|/ is multiplied by 180 and divided by n. Thus the Eq. 19, vj/ = 21.13. Since nl>0 the v|/ = -21.13. 

CONCLUSIONS 


The dihedral angle ‘v|/’ of an amino acid is calculated, which can be the input for rotating the dihedral angle with 
respect to the bond length of atoms T, rotation angle ‘0’ and the dihedral angle ‘vj /’. These 3 values can be used for 
rotation of dihedral angle to predict the structure of a protein. 
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